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ABSTRACT 


f 


A  difference  method*  which  for  one  choice  of  weights  applied  to  the 
differences  leads  to  -qually  weighted  least  squares  polynomial  fits 
to  equally  spaced  data,  is  generalized  so  that  the  approximation  can 
he  obtained  in  terms  of  polynomials  chosen  by  the  user  in  the  fitting 
of  equally  spaced  or  unequally  spaced  data. 

FOREWORD 


Using  tills  method,  an  investigation  was  made  into  the  ballistics 
field  with  particular  emphasis  on  the  preparation  of  bombing  tables. 
The  work  was  conducted  under  Task  Assignment  Number  RMM 0-42007/2101/ 
F00809001.  Credit  is  due  Richard  Hageman,  Richard  Shannon,  Gordon 
Barker,  and  Mi's.  Jane  H.  Martina  for  their  assistance  in  the  prepara¬ 
tion  of  the  tables  found  in  Appendix  A  and  the  example  in  Appendix  B. 

INTRODUCTION 

The  technique  of  using  differences  in  making  polynomial  fits  to 
discrete  data  is  not  new  and  has  been  widely  discussed  in  the  litera¬ 
ture.  However,  the  method  described  in  this  report  enables  the  user 
to  apply  this  method  using  any  form  of  interpolating  polynomial  which 
best  suits  his  purpose.  Also,  by  changing  from  ordinary  differences, 
which  are  used  in  working  with  equally  spaced  data,  to  divided  differ¬ 
ences,  the  method  can  be  used  for  unequally  spaced  data. 


For  the  convenience  of  the  user  Appendix  A  contains  the  tabulations  of 
the  multipliers,  the  weights  and  the  summations  of  weights  through  the 
4th  difference  for  equally  spaced  data  containing  from  5  to  30  points. 


With  the  modification  presented  in  the  text  on  Page 
may  be  adapted  for  use  with  unequally  spaced  data. 


9  these  weights 


EQUALLY  SPACED  DATA 

Exact  nth  Degree  Fit  to  n  +  1  Points 

Newton,  Gauss,  and  others  have  given  formulas  for  obtaining  an  nth 
degree  fit  to  n  +  1  equally  spaced  points  (sec  reference  (a)).  Such 
formulas  have  the  advantage  of  being  readily  available,  but  the  dis¬ 
advantage  of  not  giving  the  result  in  the  form  desired  for  many  appli¬ 
cations.  At  the  expense  of  additional  computation  it  is  possible  to 
get  the  interpolating  polynomial  in  the  desired  form  for  any  applica¬ 
tion  . 


1 


Suppose  this  form  for  iho  nth  degree  polynomial  through  n  f  1  points 
is 


fn  =  A0P0  +  Ai?!  +  A2P5  +  . . .  +  AnPn  .  (1) 

P]c  can  he  any  kth  degree  polynomial.  Tire  argument  used  in  P,c  can  be 
the  original  independent  variable,  x ,  or  some  other  argument,  say  u, 
which  is  a  linear  function  of  x, 

Tf  the  n  t-  1  points  are  substituted  in  equation  (1).  n  +  1  equations 
in  A0,  Alf  ...  Ar  are  obtained.  Since 

AnPj  =  0  for  J  <  n  ,  (2) 

A“f  n  -  AnAnPn 

However, 

ArPr  «  h**r  <ri  (3) 

where  Pr  <r’  is  the  rth  derivative  of  Pr  with  respec*-  to  the  argument 
used,  and  h  is  the  increment  in  that  argument  between  successive 
points.  Equation  (2)  can  then  be  solved  for  Au,  with 


Then 


=  A"fn  -  &nfn 
”  AnPn  ’  hl¥n  ,n> 


f 


n-i 


=  f. 


Anfa  p 
hnPn  41  > 


(4) 


will  be  a  polynomial  of  degree  n  -  !• 

If  the  first  n  values  of  fn_3  ore  used  it  is  now  possible  to  find 


An-i 


j^n-lp  tn-1) 
n-i 


(5) 


and 


An-if  . 

f  =  f  ~  - -  P  ,  (6) 

n-2  n-1  hn~lp  'Ji-1’  n_1  . 

n-i 


2 


will  be  a  polynomial  of  degree  n  -  2 ,  The  first  n  -  1  vainer.  !' 

f  „  can  be  used  to  find  A  and  so  on  until  all  the  A’s  have  beer, 
n-a  n-;j 

evaluated .  The  Newton  formula  uses  P  'a  chosen  in  ouch  a  way  that 

Pn  vanishes  for  the  n  points  used  in  finding  A..  ,,  P  und  I 

vanish  for  the  i:  -  1  points  used  in  finding  A„ etc.,  which  is 
ideal  provided  the  form  obtained  for  f  »  is  convenient  for  use.  For 

some  purposes  it  is  desirable  to  have  Pk  -  x.  It  is  obvious  in 

such  a  case  thut  the  derivatives  and  integrals  may  be  seen  at  a 

glance  if  necessary. 

Another  purpose  might  call  for  the  use  of  the  Tchebyehev  polynomials 
on  the  right  hand  side  of  equation  (.1)  if  truncation  to  a  lower  degre 
appeal’s  desirable.  Other  uses  of  fn  will  demand  other  forms  of  Pk, 

Smoothed  nth  Degree  Fit  to  n  +  1  +  m  Points 

Let  it  be  assumed  that,  the  function  defined  by  the  n  +  1  +  m  values 
is 


r  “  f„  +  e  (7) 

where  f n  is  an  nth  degree  polynomial  and  e  is  the  error  of  ob servo' i > 

or  noise.  If  a  satisfactory  method  for  computing  (Aky),  the  mean  kth 

difference,  is  available  and  it  can  be  assumed  that  (Ake)  will  be 
sufficiently  small  for  all  k,  then  the  approximation  of  r, 

f  n  =  AoP0  +  A1P1  +  ASPS  +  *  4  *  +  AnPn  > 


can  be  found. 


Since  c  is  involved  it  will  be  necessary  to  use  n  +  1  +  m  values  in 
the  computation  of  every  mean  difference.  An  can  be  found  from 


Let 


/  V 

(A“f)  AnA"pn  ,An^hnPn'n,j  . 


A 

f„ 

.  (• 

r  *  e 
n 

A 

f  , 

-  'r 

A  P  -  i*  ,  »  e 

n-l 

n 

n  n  n-l 

A 

f 

/> 

-  1  « 

> 

ii 

n-r* 

n-l 

n-i  n-i  n-2 

t  e 


etc , 


(0) 
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Then,  for  ouch  Ak, 


compcu'1  son  of  equations  (11)  and  (12)  allows  that.  Uh*  m'a  are  minus  the 
first  differences  of  the  act  (Ow^WgWjj...  w  0) .  Given  the  v’b,  the  m  *  r. 
con  be  found  and  given  the  m'e,  the  w's  can  be  found.  All  the 
weighted  mean  first  differences  co*i  be  found  by  these  formulas. 


There  are  many  ways  in  which  the  m'e  (or  w'a)  oon  be  chosen.  One 
method  of  molting  the  choice  is  to  minimize  the  maximum  absolute  value 

111 1 

of  —  (see  reference  (b)).  A  choice  which  is  more  consistent  with 

practice  is  to  assume  the  individual  err ora  to  be  uncorre luted  and  to 
have  equal  variance.  With  this  in  mind  let  the  w's  be  chosen  in  such 

a  way  as  to  minimize  M  =  £  m2>  holding  W  =  £w  constant,  w.hich  will 
minimize  the  variance  of  ( Zy)  . 


Then 


(1M  -  2 


m. 


mr 


dmr 

dw£ 


dm<> 

+  llli 

dw^ 

dw^ 

dn^ 

+  m. 

-irTJ1 

+  m. 

X 

dw  a 

dw. 


“a 

(JWg 


dw. 


(13  ) 


dmr.0 

4-  m  J  v. 

^  ^»r-i  ,t 

dm ,, 

r»ij,  .  _ 

1  .  1 
j  ciw,,j 

dw 

r 

dw 

1' 

dwr 

dW  -  dw1  +  dws  -t  dw  r  »  0 


a  <t  > 


The  requirement  for  a  constrained  minimum  is  met  if  the  rank  of  the 
matrix  of  the  coefficients  of  the  dw's  in  equations  (13)  and  (it)  is 
one.  Therefore,  the  coefficients  of  the  dw's  In  equation  (13)  are  all 
equal.  From  equations  (11)  and  (12)  each  of  these  coefficients  is  a 
constant  times  a  first  difference  of  the  m's.  The  in's,  in  turn,  are 
mir.ue  the  first  differences  of  the  set 


(Ow  wn . , ,wr0] 

o  that  this  seL  has  a  constant  second  difference  and  the  w  function 
:s  a  quadratic  with  two  known  zeros. 


S 


The  weigh  leu  mean  of  the  second  r{  i  (Terence  can  be  writ1. .ess 


or 


1 
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0 
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, 

•  »  ♦  » 

• 
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* 

• 

• 

> 
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• 

* 

0 

0 

0 
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0 
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yr 

0 

0 

0 

0 

0 

.  .0  .1 

n 

—  c. 

0 

0 

0 

0 

0 

o 

q 

1 

(A3y)  •-  (m^y  o  ■' 

mxyx 

msy2 

«  « 

m  i\V  r) 

(is) 

’Hie  m'a  are  the  second  differences  of  tlv»  set 

[  OOw  1w,, .  .  w  (,_  jOO] 


constant  i&uuu  to  the  roqui  rement 


and  minimizing  S  with  S  w  held  c 

that  the  w  Arne  lion  be  a  fourth  degree  polynomial  with  four  Known 
zeros  * 


In  finding  the  mean  klh  difference,  (/l'y)>  the  »''s  will  be  («l)k 
times  the  It  th  difference  of  the  set 

lOO. . .0w1w2...wr_k+10...00) 

having  k  zeros  at  each  end.  The  requirement  that  £  m‘?  be  a  minimum 
with  £  w  held  constant  leads  to  a  2k th  degree  polynomial  with  2k  known 
zeros  for  the  w  function. 
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The  ith  weight  in  getting  the  mean  kth  difference  from  r  +  1 
points  can  thus  be  written 


where  p  is  some  arbitrary  constant.  The  weights  given  in  Table  2  in 
Appendix  A  were  obtained  by  using  equation  (16).  It  is  seen  from  the 
equation  that  for  k  =  0,  p  Wj_  r  1  for  all  i. 

Since  the  m's  for  finding  the  mean  kth  difference  satisfy  a  kth  degree 
polynomial  and  since  the  (k  +  j  )th  difference  of  a  kth  degree  poly¬ 
nomial  is  zero  for  j  >  0,  the  m  polynomials  are  orthogonal  over  an 
equally  spaced  set  with  equal  weighting  and  are  identical,  except 
for  an  arbitrary  multiplier,  with  orthogonal  polynomials  for  equal 
spacing  and  equal  weighting  discussed  on  pages  287  -  291  of  refer¬ 
ence  (a)  and  given  on  pages  375  -  381  of  reference  (c).  The  use  of 
equation  (16)  in  finding  the  mean  kth  difference  thus  leads  to  an 
equally  weighted  least  squares  polynomial  fit  of  equally  spaced  data 
for  all  possible  choices  of  the  IV s  in  equation  (1). 

The  results  obtained  using  the  w  polynomials  of  equation  (IS)  and  the 
resulting  in  polynomials  for  the  P^'s  in  equation  (l)  are,  of  course, 
very  familiar.  The  use  of  the  w's  to  fit  in  terms  of  general  P^'s  and 
the  modification  of  the  w's  for  fitting  unequally  spaced  data  are  not 
so  familiar.  More  general  methods  for  choosing  the  w’s  and  m's  may 
well  be  a  rewarding  field  for  further  study. 

UMEQPALLY  SPACED  DATA 

Exact  nth  Degree  Fit  to  n  +  1  Points 


The  method  for  equally  spaced  data  may  be  modified  for  use  with 
unequally  spaced  data  by  using  divided  differences  instead  of  ordinary 
differences . 

Thus,  if 


f  =  A  P  +  A,F,  +  A  P  +  ...  +  A  P  , 
n  oo  11  22  n  n  7 


■K,  V  Un]  =AnPntuo>  V 


un3 


(i) 

(17) 
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where 


f[V  u 


i> 


= 


f(u0) 


(u0  - 


ui> 


(u0  -  un) 


f(ui) 

(ux  -  u0)  .  .  ■  (ux  -  un) 


f(un) 


(Un  -  u0>  ...  (Un  -  un.x) 

for  any  positive  integer  n,  is  the  nth  divided  difference  defined  in 
reference  (a). 

If 

Pn(u)  =  B0  +  Bxu  +  B2u2  +  . . .  +  Bnun  , 

the  nth  divided  difference  of  Pn(u)  is 

P  [  u  ,  vi , ,  u_  . . .  u  ]  =  B 
n  o'  i*  2  n  n 


So  that,  except  for  the  arithmetic  required  in  the  compute. tiun  of 
ftu0,  Uj_,  ug,  ...  u  ]  ,  An  is  easily  found  from  equation  (17). 

The  finding  of  f  ,  f  ,  — ,  and  A  ,  A  ....  follows  a  path 
°  n-j/  1  n-2'  '  n-.l'  n-2 

parallel  to  thao  described  for  equally  spaced  data. 

Smoothed  nth  Degree  Fit  to  n  +  1  +  m  Points 


Let 


(W'f) 


be  the  mean  kth  divided  difference  of  f.  Then  equation  (9)  can  be 
rewritten 


Ak 


(18) 


and  the  problem  of  fitting  is  then  solved  if  a  satisfactory  way  can 
be  found  for  finding  the  mean  kth  divided  difference. 
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Finding  Mean  Divided  Differences 


The  mean  kth  divided  difference  can  be  found  from 


<Akf) 


(.0 


i  ftu0,  ul> 


uk]  +  U>2  ftUjL,  U2,  ...  Uk+1J  + 


.  (19) 


The  optimum  choice  of  the  oi's  is  a  more  difficult  problem  than  was  the 
choice  of  the  w's  in  the  case  of  equally  spaced  data.  A  possible 
choice  is  to  take 


ui  =  <uk  -  uo>  wi 
=  (uic+i  -  ux) 


wl  =  (uk+i+x  ”  ui+x)  wi  > 


the  w's  being  those  chosen  for  finding  the  mean  kth  difference  for 
equally  spaced  data.  For  finding  the  differences  when  k  -  0,  choose 
w1  =  i , 

This  will  approach  the  cquallv  spaced  fitting  as  the  intervals 
become  more  nearly  equal  and  may  be  expected  to  give  reasonable 
results  in  other  cases. 

It  should  be  mentioned  that  the  interesting  difference  relation 
between  the  w's  and  m's,  and  the  resulting  generation  of  the  m  poly¬ 
nomials,  for  equally  spaced  data  does  not  hold  true  for  the  w's  and 
m's  used  with  unequally  spaced  data  involving  divided  differences . 

Given  in  Appendix  B  is  an  example  using  this  divided  difference 
technique  in  smoothing  y  -  sin  x  over  a  non-constant  data  interval. 
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APPENDIX  A 


TABLE  1 


MULTIPLIERS  FOR  EQUALLY  SPACED  DATA 


5  POINTS 


i 

m(0) 

m(l) 

m<2) 

m(3) 

m(4) 

0 

1 

*■2 

2 

-1 

1 

1 

1 

••1 

-1 

2 

-4 

2 

1 

0 

-2 

0 

6 

3 

1 

1 

•I 

-2 

-4 

4 

1 

2 

2 

1 

1 

Zw 

5 

10 

7 

2 

1 

E<m)2 

5 

10 

14 

10 

70 

6  POINTS 

i 

m(0) 

m(l> 

m(2) 

m(3) 

m(4) 

0 

1 

2 

3 

4 

5 


1 

1 

1 

1 

1 

1 


-5 

-3 

-1 

1 

3 

5 


5 

_  ^ 

■“  a. 

-4 

-4 

*1 

5 


-5 

7 

4 
-4 

-7 

5 


1 

-3 

2 

2 

i 


£v 

£<m)2 


6 

6 


35 

70 


28 

84 


18 

180 


2 

28 


1 


TABLE  1  (Continued) 


Q  POINTS 


i 

m(O) 

m(l) 

m(2> 

m  (3) 

m(4) 

0 

1 

-7 

7 

-7 

7 

1 

1 

-L 

-5 

i 

5 

=  13 

2 

1 

-3 

•>3 

7 

-3 

T 

V 

1 

-1 

-5 

3 

9 

4 

1 

1 

-5 

-3 

9 

5 

1 

3 

-3 

-7 

-3 

6 

1 

5 

1 

-5 

-13 

7 

1 

7 

7 

7 

7 

Q 

84 

04 

66 

44 

£«a 

8 

168 

168 

264 

616 

2 


TABLE  1  (Continued  ) 
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TABLE’  1  (Continued) 


13  POINTS 


i 

m(0) 

m(l) 

m  (a) 

m(3) 

m(4) 

0 

1 

-6 

22 

-11 

99 

1 

1 

*>5 

11 

0 

-66 

2 

1 

-4 

2 

6 

-96 

3 

1 

-3 

“5 

Q 

-54 

4 

1 

-2 

-10 

7 

11 

5 

1 

-1 

-13 

4 

64 

6 

1 

0 

-14 

0 

84 

7 

1 

1 

-13 

-4 

64 

8 

1 

2 

-10 

-7 

11 

9 

1 

3 

-5 

-e 

-54 

10 

1 

4 

2 

-6 

-96 

11 

1 

5 

11 

0 

-66 

12 

1 

6 

22 

11 

99 

71  w 

13 

182 

1001 

572 

4862 

13 

182 

2002 

572 

68,068 

5 


TABLE  1  (Continued) 


14  POINT'S 


i 

i»  00 

m(l) 

«<2) 

•  (  ) 

m(4) 

0 

1 

-13 

13 

-  .i. :  3 

143 

1 

1 

-11 

7 

-]  1 

-  77 

o 

c. 

1 

-9 

2 

uo 

-132 

3 

1 

-7 

-2 

:.'C 

-92 

4 

1 

-5 

-5 

95 

-13 

5 

1 

-3 

-7 

C7 

63 

6 

1 

-1 

-8 

'  1  4 
rJ 

100 

7 

1 

1 

-0 

-24 

108 

0 

1  1 

3 

-7 

-G7 

63 

9  I 

!  i 

5 

-5 

-35 

-13 

10 

1  i 

7 

-2 

-98 

-92 

11 

:  i 

9 

2 

-66 

-132 

1” 

i 

11 

7 

11 

-77 

13 

i 

13 

13 

143 

143 

i 

14 
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728 

.724 

9724 

Z<  )'  1 

14 
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728 

9724 

136,136 
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m(0) 

in  (1  > 

m<2) 

m<3) 

m  (4) 

0 

1 

-7 

91 

-91 

1001 

1 

1 

*•6 

52 

-13 

-429 

2 

1 

-5 

19 

35 

-869 

3 

1 

-4 

-8 

58 

-704 

4 

1 

»3 

-29 

61 

-249 

5 

1 

-2 

-44 

49 

251 

6 

1 

-I 

-53 

27 

621 

7 

1 

0 

-56 

0 
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6 

1 

1 

-53 

-27 

621 

9 

1 

2 

-44 

-49 

251 

10 

1 

3 

-29 

-61 

-249 

11 

1 

4 

-8 

-58 

-704 

12 

1 

5 

19 

-35 

-669 

13 

1 

6 

52 

13 

-429 

14 

1 

7 

91 

91 

1001 

15 

2  GO 

6108 

7956 

92, 378 

Z 

IS 

280 

37,126 

39, 780 

6,466,460 
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i 

m(0) 

0 
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1 

2 

1 

3 

1 

4 

1 

5 

1 

6  1 

1 

7  '  | 

1 

8 

1 

9 

1 

10 

1 

11 

1 

12 

1 

13 

1 

14 

1 

15 

1 

3.6 

1 

Z,v 

17 

17 

m(i) 

m(2) 

m(3) 

-8 

40 

-26 

-7 

25 

-7 

•*6 

12 

7 

-5 

1 

IS 

-4 

-6 

18 

-3 

-15 

17 

-2 

-20 

13 

-1 

-23 

7 

0 

-24 

0 

1 

-23 

-7 

2 

-20 

-13 

3 

-15 

-17 

4 

-8 

-18 

6 

1 

-15 

6 

12 

-7 

7 

25 

7 

8 

40 

28 

408 

3876 

3876 

408 

7752 

3876 

m(4) 


52 

-13 

-39 

-39 

-24 

-3 

17 

31 

36 

31 

17 

-3 

-24 

-39 

-39 

-13 

52 


8398 

16,796 
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20  POINTS 


i 

m(0) 

m(l) 

m<2) 

rt  (3) 

m(4) 

0 

1 

-19 

57 

-969 

1936 

1 

1 

-17 

39 

-357 

-102 

2 

1 

-15 

23 

85 

-1122 

3 

1 

-13 

9 

377 

-1402 

4 

1 

-11 

-3 

539 

-1187 

5 

1 

-9 

-13 

591 

-687 

6 

1 

-7 

-21 

553 

-77 

7 

1 

-5 

-27 

445 

503 

8 

1 

-3 

-31 

.  287 

948 

9 

1 

-1 

-33 

99 

1168 

10 

1 

1 

-S3 

-99 

1188 

11 

1 

3 

-31 

-267 

948 

12 

1 

S 

-27  • 

-445 

503 

13 

1 

7 

-21 

-553 

-77 

14 

1 

9 

-13 

-591 

-687 

15 

1 

u 

-3 

-539 

-1187 

16 

1 

13 

9 

-377 

-1402 

17 

1 

15 

23 

-85 

-1122 

18  • 

1 

17 

39 

357 

-102 

IS 

i 

*i  ^ 

5? 

969 

1238 

£w 

20 

1330 

8778  • 

245,157 

653, 752 

ZOn)2 

20 

2660 

17,656 

4,903,140 

22,881,320 

TABLE  1  (Continued) 


21  POINTS 


i 

m(0) 

m(l) 

»fe> 

m(3) 

mOi) 

A 

u 

1 

-10 

190 

-285 

969 

1 

1 

«•  9 

133 

-114 

0 

2 

1 

-8 

82 

12 

-510 

3 

1 

-7 

37 

98 

-680 

4 

1 

-6 

-2 

149 

-615 

5 

1 

-5 

-36 

170 

-406 

6 

1 

'  -4 

.62 

166 

-130 

7 

1 

-3 

•83 

142 

150 

8 

1 

-2 

•  98 

103 

385 

9 

1 

-1 

-107 

54 

540 

10 

1 

0 

-110 

0 

594 

11 

1 

1 

*•107 

-54 

540 

12 

1 

2 

•98 

-103 

365 

13 

1 

3 

•88 

-142 

150 

14 

1 

4 

•62 

-166 

-130 

15 

1 

5 

•36 

-170 

-406 

16 

x 

6 

-2 

-149 

-615 

17 

1 

7 

37 

-98 

-680 

18 

1 

8 

82 

-12 

-510 

19 

1 

q 

i  ** 

AW 

1  T  A 
11*1 

A 

20 

1 

10 

190 

285 

969 

Zv 

21 

770 

33,649 

86,526 

408,595 

2>)a 

21 

770 

201,894 

432,630 

5,720,330 
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i 

m(0) 

m(l) 

m<2) 

m(3) 

m(4) 

0 

1 

-21 

35 

-133 

1197 

1 

1 

-19 

25 

-57 

57 

2 

1 

-17 

16 

0 

-570 

3 

1 

-15 

8 

40 

-810 

4 

1 

-13 

1 

65 

-775 

6 

1 

-11 

-5 

77 

-563 

6 

1 

-9 

-10 

78 

-258 

7 

1 

-7 

-14 

70 

70 

8 

1 

-5 

-17 

55 

365 

9 

1 

-3 

-19 

35 

585 

10 

1 

-1 

-20 

12 

702 

11 

1 

1 

-20 

-12 

702 

12 

1 

3 

-19 

-35 

585 

13 

1 

5 

-17 

-55 

365 

14 

1 

7 

-14 

-70 

70 

15 

1 

9 

-10 

-78 

-258 

16 

1 

11 

-5 

-77 

-563 

17 

1 

13 

1 

-65 

-775 

10  ' 

15 

0 

-40 

-810 

19 

**  n 

X  i 

16 

0 

=570 

20 

1 

19 

25 

57 

57 

21 

1 

21 

35 

133 

1197 

2  w 

22 

1771 

7084 

48,070 

624,910 

J>>2 

22 

3542 

7084 

96,140 

8, 748, 740 
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i 

m(0) 

m  (l) 

in  (?.) 

m  (3  ) 

m  (-1 ) 

0 

1 

-12 

92 

-506 

1510 

1 

1 

-11 

69 

-253 

253 

2 

1 

-10 

40 

-55 

-517 

3 

1 

-9 

29 

93 

-897 

4 

1 

-0 

12 

196 

5 

1 

-7 

*3 

259 

-857 

6 

1 

-6 

-16 

287 

-597 

7 

1 

*•5 

-27 

285 

-267 

8 

1 

-4 

*36 

258 

78 

9 

1 

-3 

-43 

211 

393 

10 

1 

-2 

-48 

149 

643 

11 

1 

-1 

-61 

77 

003 

12 

.1 

0 

-52 

0 

858 

13 

1 

1 

-51 

-77 

803 

14 

1 

2 

-48 

-149 

643 

15 

1 

3 

-43 

-211 

393 

16 

1 

4 

-36 

-250 

78 

17 

1 

6 

-27 

-285 

-267 

10 

1 

6 

-16 

-287 

-597 

19 

1 

7 

-3 

-259 

-857 

20 

1 

0 

12 

-196 

-  982 

21 

1 

9 

29 

-93 

-097 

22 

1 

10 

48 

55 

-517 

23 

1 

11 

69 

253 

253 

24 

1 

12 

02 

506 

1518 

E* 

25 

1300 

26, 010 

296,010 

1,430,715 

E(m)2 

25 

1300 

53,820 

1,480,050 

14,307,150 
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i 

m(0) 

m(l) 

m(&) 

m(5) 

m  (4 ) 

0 

1 

-2G 

50 

-3450 

2530 

1 

1 

-23 

38 

-1794 

506 

p 

1 

*•21 

27 

*403 

-759 

3 

1 

-19 

17 

513 

-1419 

4 

1 

-1? 

0 

1224 

“1614 

5 

1 

*15 

0 

1680 

“1470 

6 

1 

-13 

-7 

1911 

*1099 

7 

1 

*11 

*15 

1947 

-599 

8 

1 

-9 

•  10 

1010 

*54 

0 
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-7 

*22 

1664 

466 

10 

1 

-5 

*25 

1185 

905 

1.1 

1 

*3 

*27 

741 

1221 

12 

1 

-1 

*20 

252 

1386 

13 

1 

1 

*28 

-262 

1386 

14 

1 

3 

*27 

-741 

1221 

15 

1 

6 

*26 

*1185 

005 

16 

1 

7 

*22 

-1554 

466 

17 

1 

0 

*10 

*1818 

-54 

IQ 

1 

11 

*13 

-1947 

-599 

19 

1 

13 

*7 

-1911 

-1099 

20 

1 

15 

0 

*1680 

-1470 

21 

.1 

17 

0 

*1224 

*1614 

22 

1 

19 

17 

*515 

-1419 

23 

1 

21 

27 

483 

*759 

24 

1 

23 

38 

1794 

506 

25 

1 

25 

50 

3450 

2530 

Sw 

26 

2925 

16,360 

2,341,170 

2, 861/430 

£<m>2 

26 

5050 

16  ,  380 

70,235,100 

40,060,020 
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i 

m(0) 

m(l) 

0 

1 

*13 

1 

1 
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2 

1 
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*7 

1 
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4 

1 

*9 

5 

1 

-8 

6 

1 

-7 

7 

1 

**b 

8 

1 

-5 

9 

1 

-4 

10 

1 

-3 

11 

1 

-2 

12 

1 

*1 

13 

1 

0 

14 

1 

1 

15 

1 

2 

16 

1 

3 

1? 

1 

4 

18 

1 

5 

19 

1 

6 

20 

1 

7 

oi 

t—  a.  i 

1 

ml. 

O 

t_> 

22 

1 

9 

23 

1 

10 

24 

1 

11 

25 

I 

12 

26 

1 

13 

27 

1630 

I>>? 

27 

1638 

m<2> 

m(3) 

m  (4) 

325 

-130 

2990 

250 

-70 

690 

181 

*22 

*782 

no 

15 

*1587 

61 

42 

-1872 

10 

60 

-1770 

*35 

70 

-1400 

-74 

73 

-067 

-107 

70 

-262 

-134 

62 

333 

-155 

50 

070 

*170 

35 

.•285 

-179 

18 

1548 

-182 

0 

1630 

-179 

*10 

1540 

-170 

-35 

1285 

-155 

*50 

070 

-134 

-62 

330 

-107 

-70 

*262 

-74 

-73 

-067 

-35 

-70 

-1400 

» r> 

t.KJ 

/*/'■* 

**ow 

•.1  f  <U 

61 

•»  42 

- 1072 

110 

*15 

-1507 

181 

22 

-7G2 

250 

70 

690 

325 

13C 

2990 

118,755 

101, 790 

4,032,015 

712,530 

101, 790 

56,440,210 
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i 

m(0) 

in  (l) 

m(2) 

m  (3) 

m  (4  ) 

0 

1 

-27 

117 

-1755 

.1 755 

1 

1 

-25 

91 

-975 

455 

2 

1 

-23 

67 

-345 

-395 

3 

1 

-21 

45 

147 

-879 

4 

1 

“19 

25 

513 

-1074 

5 

1 

-17 

7 

765 

-1050 

6 

1 

-15 

-9 

015 

-070 

7 

1 

-13 

-23 

975 

-590 

8 

1 

-11 

-35 

957 

-259 

9 

1 

-9 

-45 

873 

81 

10 

1 

-7 

-53 

735 

395 

11 

1 

-5 

-59 

555 

655 

.12 

1 

-3 

-63 

345 

84C 

13 

1 

-1 

-65 

117 

936 

14 

1 

1 

-65 

-117 

936 

15 

1 

3 

-63 

•>  345 

840 

16 

1 

5 

-59 

-555 

655 

17 

l 

7 

-53 

-735 

395 

18 

1 

9 

-45 

-873 

81 

19 

1 

11 

-35 

-957 

-259 

20 

1 

13 

-23 

-975 

-590 

21 

1 

15 

-9 

-915 

-870 

22 

1 

17 

7 

-765 

-1050 

23 

1 

19 

25 

-513 

-1074 

24 

•4. 

21 

45 

-147 

-079 

25 

1 

23 

67 

345 

-393 

26 

l 

25 

91 

975 

455 

27 

1 

27 

11' 

1755 

1755 

28 

3G54 

47,502 

1,577,745 

? , 804 , 890 

Z>)2 

20 

7300 

05,004 

10,5<.:,U40 

19,034,160 
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i 

m(0) 

m(l) 

m  (2  ) 

m  (3) 

n  ■  ' 

0 

1 

-14 

126 

-819 

4095 

1 

1 

-13 

99 

-468 

1170 

2 

1 

-12 

74 

-182 

-780 

3 

1 

-11 

51 

44 

-1930 

4 

1 

-10 

30 

215 

-2441 

5 

1 

-9 

11 

336 

-2460 

6 

1 

-8 

-6 

412 

-2120 

7 

1 

-7 

-21 

448 

-1540 

8 

1 

-6 

-34 

449 

-825 

9 

1 

-5 

-43 

420 

-66 

10 

1 

-4 

-54 

366 

660 

11 

1 

-3 

-61 

292 

3290 

12 

1 

-66 

203 

1775 

13 

1 

-i 

-69 

104 

2080 

14 

1 

0 

-70 

0 

2184 

15 

1 

1 

-69 

-104 

2080 

16 

1 

2 

-66 

-203 

1775 

17 

n 

J. 

3 

-61 

-292 

1290 

18 

1 

4 

-54 

-366 

660 

19 

1 

5 

-45 

-420 

-66 

20 

1 

6 

-34 

-449 

-825 

21 

1 

7 

-21 

-448 

-1540 

22 

n 

-L 

8 

-6 

-4.12 

-2120 

23 

1 

9 

11 

-336 

-2460 

24 

1 

10 

30 

-215 

-2441 

25 

1 

11 

51 

-44 

-1930 

26 

1 

12 

74 

182 

-780 

27 

-L 

13 

99 

468 

1170 

28 

1 

14 

126 

819 

4095 

£w 

29 

2030 

56,637 

841,464  7 

,713, 420 

XOn)* 

29 

2030 

113,274  4, 

.207,320  107 

,987,880 
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4 

m(0) 

m(l) 

m<2) 

m(3) 

m(4) 

0 

I 

-29 

203 

-1827 

23751 

1 

1 

-27 

161 

-1071 

7371 

2 

1 

-25 

122 

-450 

-3744 

3 

1 

-23 

86 

46 

-10504 

4 

1 

-21 

53 

427 

-13749 

5 

1 

-19 

23 

703 

-14249 

6 

1 

-17 

-4 

884 

-12704 

7 

1 

-15 

-28 

980 

-9744 

8 

X 

-13 

-49 

1001 

-5929 

9 

1 

-11 

-67 

957 

-1749 

10 

1 

-82 

858 

2376 

11 

1 

-7 

-94 

714 

6069 

12 

1 

-5 

-103 

535 

9131 

13 

1 

-3 

-109 

331 

11271 

14 

1 

-1 

-112 

112 

12376 

15 

1 

1 

-112 

-112 

12376 

16 

1 

3 

-109 

-331 

11271 

-»  *7 
-J-  ( 

1 

5 

-103 

-535 

9131 

18 

i 

7 

-94 

-714 

6069 

19 

1 

9 

-82 

-858 

2376 

20 

1 

11 

-67 

-957 

-1749 

21 

1 

13 

-49 

-1001 

-5929 

22 

1 

15 

-28 

-980 

-9744 

23 

1 

17 

-4 

-884 

-12704 

24 

1 

19 

23 

-703 

-14249 

25 

1 

21 

53 

-427 

-13749 

26 

1 

23 

86 

-46 

-10504 

27 

1 

25 

122 

450 

—  3  7  4  4 

28 

1 

27 

161 

1071 

7371 

29 

1 

29 

203 

1327 

23751 

Z  V 

i 

30 

4495 

100, 688 

2,136,024 

52,451,256 

Io«>2 

!  30 

Q990 

302,064 

21,360,240 

3,671,587,920 
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5  POINTS  6  POINTS 


7  POINTS 

8  POINTS 

k=l 

k=2  k=3 

k  =4 

k=l 

k--£ 

k=3 

k=4 

i  pwt 

pvL  pw. 

pv<_ 

rVW, 

•  i. 

PW  j 

pv  | 

1  3 

2  5 

3  6 

4  6 

5  5 

6  3 

7 


2  v  28 


5 

1 

3 

7 

1 

i 

7 

7 

1 

10 

2 

5 

12 

15 

16 

15 

2 

12 

2 

3 

15 

20 

20 

15 

•z 

10 

1 

16 

20 

16 

7 

4 

5 

15 

15 

7 

5 

12 

7 

6 

7 

7 

8 

42 

6 

11 

84 

84 

65 

44 

h  cu  (O 
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9  POINTS  10  POINTS 


i 

k=i 

pwt 

k-2 

pwt 

k  =3 
pwx 

k=4 

pwi 

k  =1 
pwt 

k=2 

pWi 

k  =3 

pWi 

k  =4 

pWi 

i 

i 

4 

28 

14 

14 

9 

6 

42 

18 

1 

2 

7 

63 

35 

35 

16 

14 

112 

50 

2 

3 

9 

90 

50 

45 

21 

21 

175 

75 

3 

4 

10 

100 

50 

35 

24 

25 

200 

75 

4 

5 

10 

90 

35 

14 

25 

25 

175 

50 

5 

6 

9 

63 

14 

24 

21 

112 

18 

6 

7 

7 

28 

21 

14 

42 

7 

8 

4 

16 

6 

8 

9 

9 

9 

10 

60 

462 

198 

143 

165 

132 

858 

286 

11  POINTS 


12  POINTS 


i 

k=T 

k=2 

pwl 

k  -3 

PW< 

k=4 

pWi 

k=i 

P*i 

k  =2 
pwA 

k=3 

Pwi 

k  =4 

P«i 

i 

1 

5 

15 

30 

6 

11 

55 

33 

33 

1 

2 

9 

36 

84 

18 

20 

135 

96 

105 

2 

3 

12 

56 

140 

30 

27 

216 

168 

189 

3 

4 

14 

70 

175 

35 

32 

280 

224 

245 

4 

5 

15 

75 

175 

30 

35 

315 

245 

245 

5 

6 

15 

70 

140 

18 

36 

315 

224 

189 

6 

7 

14 

56 

84 

6 

35 

280 

168 

105 

7 

8 

12 

36 

30 

32 

216 

96 

33 

8 

9 

9 

15 

27 

136 

33 

9 

10 

5 

20 

55 

10 

11 

11 

11 

12 

£v 

110 

429 

858 

143 

286 

2002 

1287 

1144 

TABLE  2  (Continued) 


13  POINTS  14  POINTS 


i 

k-1 

k=2 

pWi 

k=3 

pVi 

k=4 

pWl 

k=l 

pWl 

k  *2 

pWi 

k  =3 

Pwi 

k  =4 

Pwi 

i 

1 

6 

22 

11 

99 

13 

13 

143 

143 

1 

2 

11 

55 

33 

330 

24 

33 

440 

495 

2 

3 

15 

90 

60 

630 

55 

55 

825 

990 

3 

4 

10 

120 

84 

882 

40 

75 

1200 

1470 

4 

5 

20 

140 

98 

980 

45 

90 

1470 

1764 

5 

6 

21 

147 

98 

882 

48 

98 

1568 

1764 

6 

7 

21 

140 

84 

630 

49 

98 

1470 

1470 

7 

0 

20 

120 

60 

330 

48 

90 

1200 

990 

8 

9 

10 

90 

33 

99 

45 

75 

825 

495 

9 

10 

15 

55 

11 

40 

55 

440 

143 

10 

11 

11 

22 

33 

33 

143 

11 

12 

6 

24 

13 

12 

13  • 

13 

13 

14 

£v 

102 

1001 

572 

4862 

455 

728 

9724 

9724 

15  POINTS  16  POINTS 


i 

k“l 

pw1 

k=2 

pw1 

k=3 

pw, 

k  =4 
pwt 

k=l 

pw1 

k=2 

pw1 

k=3 

Pw! 

k  =4 

pWj 

i 

1 

7 

91 

91 

1001 

15 

35 

455 

1365 

1 

2 

13 

234 

286 

3575 

28 

91 

1456 

5005 

2 

3 

18 

396 

550 

7425 

39 

156 

2860 

10725 

3 

4 

22 

550 

825 

11550 

48 

220 

4400 

17325 

4 

5 

25 

675 

1050 

14700 

55 

275 

5775 

23100 

5 

6 

27 

756 

1176 

15876 

60 

315 

6720 

26460 

6 

7 

28 

784 

1176 

14700 

63 

336 

7056 

26460 

7 

0 

28 

756 

1050 

11550 

64 

336 

6720 

23100 

8 

9 

27 

675 

825 

7425 

63 

315 

5775 

17325 

Q 

10 

25 

550 

550 

3575 

60 

275 

4400 

10725 

10 

11 

22 

396 

286 

1001 

55 

220 

2860 

5005 

11 

12 

16 

234 

91 

48 

156 

1456 

1365 

12 

13 

13 

91 

39 

91 

455 

13 

14 

7 

28 

35 

14 

15 

15 

15 

16 

Zw 

280 

6188 

7956 

92,378 

680 

2856 

50, 386 

167,960 

TABIiE  2  (Continued) 


17  POINTS  IQ  POINTS 


i 

k=l 

PW, 

k=2 

pw, 

k  =3 
pw. 

k-4 

PW, 

k-1 

pw, 

k=2 

pw, 

k=3 

pw, 

k  =4 
pw. 

i 

i 

8 

40 

28 

52 

17 

68 

68 

68 

1 

2 

15 

105 

91 

195 

32 

180 

224 

260 

2 

3 

21 

182 

182 

429 

45 

315 

455 

585 

3 

4 

26 

260 

286 

715 

56 

455 

728 

1001 

4 

5 

30 

330 

385 

99C 

65 

585 

1001 

1430 

5 

0 

33 

385 

462 

1188 

72 

693 

1232 

1782 

6 

7 

35 

420 

504 

1260 

77 

770 

1386 

1980 

7 

8 

36 

432 

504 

1188 

80 

810 

1440 

1980 

8 

9 

36 

420 

462 

990 

81 

810 

1386 

1782 

9 

10 

35 

385 

385 

715 

80 

770 

1232 

1430 

10 

11 

33 

330 

286 

429 

77 

693 

1001 

1001 

11 

12 

30 

260 

182 

195 

72 

585 

728 

585 

12 

13  • 

26 

182 

91 

52 

65 

455 

455 

260 

13 

14 

21 

105 

28 

56 

315 

224 

68 

14 

13 

15 

40 

45 

180 

68 

15 

16 

8 

32 

68 

16 

17 

17 

17 

18 

Sw 

408 

3876 

3876 

8398 

969 

7752 

11, 628 

14,212 

TABLE  2  (Continued) 


19  POINT'S 


20  POINTS 


k  *1  k=2  k  =3 

pwt  pwt  pWj 


1 

9 

51 

204 

2 

17 

136 

680 

3 

24 

240 

1400 

4 

30 

350 

2275 

5 

35 

455 

3185 

6 

39 

546 

4004 

7 

42 

616 

4620 

8 

44 

660 

4950 

9 

45 

675 

4950 

10 

45 

660 

4620 

11 

44 

616 

4004 

12 

42 

546 

3185 

13  • 

39 

455 

2275 

14 

35 

350 

1400 

15 

30 

240 

680 

16 

24 

136 

204 

17 

17 

51 

.18 

O' 

v/ 

19 

Iw 

570 

6783 

42,636 

k=4 

pwt 

k  -1 

pWj 

k  ==2 
pw  i 

k  =3 

(J*l 

k  =4 

fJwi 

i 

612 

19 

57 

969 

1930 

1 

2300 

36 

153 

3264 

7650 

2 

5460 

51 

272 

6000 

17850 

3 

9555 

64 

400 

11200 

31850 

4 

14014 

75 

525 

15925 

47775 

5 

18018 

84 

637 

20304 

63063 

6 

20790 

91 

728 

24024 

75075 

7 

21780 

96 

792 

26400 

Q1675 

8 

20790 

99 

825 

27225 

81675 

9 

16018 

100 

825 

2 6400 

75075 

10 

14014 

99 

792 

24024 

63063 

11 

9555 

96 

728 

20384 

47775 

12 

5460 

91 

637 

15925 

31850 

13 

2380 

04 

525 

11200 

17850 

14 

612 

75 

400 

6800 

7650 

15 

64 

272 

3264 

1938 

16 

51 

153 

969 

17 

56 

57 

18 

19 

19 

20 

163, 430 

1330 

877Q 

245,157 

653, 752 

5 


TABLE  2  (Continued) 


21  POINTS  22  POINTS 


i 

k=l 

pw: 

k=2 

PWj^ 

k=3 

pw1 

k=4 

pw. 

k=l 

pvi 

k=2 

pw, 

k=3 

pwA 

k=4 

pw1 

1 

1 

10 

190 

205 

969 

21 

35 

133 

1197 

1 

2 

19 

513 

969 

3076 

40 

95 

456 

4845 

2 

3 

27 

910 

2040 

9180 

57 

171 

969 

11620 

3 

4 

34 

1360 

3400 

16660 

72 

255 

1632 

21420 

4 

5 

40 

1000 

4900 

25480 

05 

340 

2380 

33320 

5 

6 

45 

2205 

6370 

34398 

96 

420 

3136 

45864 

6 

7 

49 

2548 

7644 

42042 

105 

490 

3822 

57330 

7 

0 

52 

2808 

8580 

47190 

112 

546 

4360 

66066 

0 

9 

54 

2970 

9075 

49005 

117 

565 

4719 

70785 

9 

10 

55 

3025 

9075 

47190 

120 

605 

4040 

70785 

10 

11 

55 

2970 

8500 

42042 

121 

605 

4719 

66066 

11 

12 

54 

2808 

7644 

3439Q 

120 

565 

4360 

57330 

12 

13  . 

52 

2548 

6370 

25400 

117 

546 

3622 

45064 

13 

14 

49 

2205 

4900 

16660 

112 

490 

3136 

33320 

14 

15 

45 

1800 

3400 

9180 

105 

420 

2380 

21420 

15 

16 

40 

1360 

2040 

3876 

96 

340 

1632 

11628 

16 

17 

34 

918 

969 

969 

85 

255 

969 

4845 

17 

18 

27 

513 

285 

i  <1 

171 

456 

1197 

10 

IS 

19 

190 

57 

95 

153 

19 

20 

10 

40 

35 

20 

21 

21 

21 

22 

Zw 

770 

33, 649 

06,526 

408,595 

1771 

7084 

48,070 

624,910 

TABLE  2  (Continued) 


23  POINTS 


i 

k11! 

PWj, 

k=2 

pvt 

k  =3 

P*i 

1 

11 

77 

77 

2 

21 

210 

266 

3 

30 

380 

570 

4 

38 

570 

969 

5 

45 

765 

1428 

6 

51 

952 

1904 

7 

56 

1120 

2352 

8 

60 

1260 

2730 

9 

63 

1365 

3003 

10 

65 

1430 

3146 

11 

66 

1452 

3146 

12 

66 

1430 

3003 

13 

.  65 

1365 

2730 

14 

63 

1260 

2352 

15 

60 

1120 

1904 

16 

56 

952 

142Q 

17 

51 

765 

969 

18 

45- 

570 

570 

IS 

38 

380 

266 

20 

30 

210 

77 

21 

21 

77 

22 

11 

23 

£v 

1032 

17, 710 

32,890 

k  =4 

pVj. 

kal 

PWi 

k<2 

P*! 

1463 

23 

253 

5985 

44 

693 

14530 

63 

1260 

27132 

80 

1900 

42840 

95 

2565 

59976 

108 

3213 

76440 

119 

3808 

90090 

128 

4320 

99099 

135 

4725 

102245 

140 

5005 

99099 

143 

5148 

90090 

144 

5140 

76440 

143 

5005 

59976 

140 

4725 

42840 

135 

4320 

27132 

128 

3808 

14535 

119 

3213 

5985 

100 

2565 

1463 

95 

1900 

80 

1260 

63 

693 

44 

253 

23 

937,365 

2300 

65, 780 

24  POINTS 


It  =3 

k  =4 
pw1 

i 

1771 

253 

1 

6160 

1045 

2 

13300 

2565 

3 

22800 

4845 

4 

33915 

7752 

5 

45696 

11016 

6 

57120 

14280 

7 

67200 

17160 

8 

75075 

19305 

9 

80080 

20449 

10 

81796 

20449 

11 

80080 

19305 

12 

75075 

17160 

15 

67200 

14280 

14 

57120 

11016 

15 

45696 

7752 

16 

33915 

4845 

17 

22800 

«•  V  w 

^  a 

13300 

1045 

19 

6160 

253 

20 

1771 

21 

22 

23 

24 

880,030  197,340 


7 


'IABLE  2  (Continued ) 


2(5  POINTS  26  POINTS 


i 

k  =1 
pwx 

k=2 

PWj 

k=3 

pWi 

k=4 

pw1 

k“l 

pwt 

k  =2 

k  =3 

k=4 

pwi 

.1 

i 

12 

92 

606 

1518 

26 

100 

2300 

2530 

1 

2 

23 

253 

1771 

6325 

48 

276 

8096 

10626 

2 

3 

33 

46t 

3050 

15G75 

69 

606 

17710 

26565 

3 

4 

42 

700 

6650 

29025 

00 

770 

30800 

51205 

4 

5 

50 

950 

9075 

40450 

105 

1050  ' 

46550 

83790 

5 

6 

57 

1197 

13566 

69760 

120 

1330 

63840 

122094 

6 

7 

63 

1428 

17136 

91000 

133 

1686 

81396 

162792 

7 

0 

68 

1632 

20400 

112200 

144 

1036 

97920 

20196C 

0 

9 

72 

1800 

23100 

128700 

163 

2040 

112200 

235620 

9 

10 

75 

1926 

25025 

139426 

.160 

2200 

323200 

260260 

10 

11 

77 

2002 

26026 

143143 

166 

2310 

130130 

273273 

11 

12 

70 

2020 

26026 

139425 

166 

2366 

132496 

273273 

12 

13  . 

70 

2002 

25025 

120700 

169 

2366 

130130 

260260 

13 

14 

77 

1926 

23100 

112200 

168 

2310 

123200 

235620 

14 

15 

75 

1800 

20400 

91800 

166 

2200 

112200 

201960 

15 

16 

72 

1632 

17136 

69768 

160 

2040 

97920 

162792 

16 

17 

60 

1420 

13566 

40450 

153 

1836 

81396 

122094 

17 

in 

63 

1197 

9975 

29925 

.144 

1596 

63840 

03790 

i  n 

-i.'w' 

19 

57 

950 

6660 

15675 

133 

1330 

46550 

51205 

19 

20 

50 

700 

3050 

6325 

320 

1050 

30000 

26565 

20 

21 

42 

462 

1771 

1518 

105 

770 

1771.0 

10626 

21 

22 

33 

253 

506 

80 

BOB 

8096 

2530 

22 

23 

23 

92 

69 

276 

2300 

23 

24 

12 

40 

loo 

24 

25 

26 

25 

26 

2 w 

1300  26,910 

l'j.3,0051 

1,430, 715 

2925 

16 >380 

1,560,700s 

1,430,715s 

AIn  computing  m(3)  in  Table  1  the  weights  used  where  2  times  these 
weights. 

2In  computing  m (3 }  in  Table  1  the  weighto  used  where  l-l/2  timee  these 
weights . 

3In  computing  m (4)  in  Table  1  the  weighto  used  where  2  times  these 
weights. 
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TABI.B  2  (Continued) 


27  POINTS  80  POINTS 


k-1 

PWi 

k"2 

pwt 

k=3 

pw1 

k=4 

pwx 

k*l 

ow1 

k~2 

PVj 

k*3 

pwt 

k-4 

pwt 

i 

13 

325 

130 

2990 

27 

117 

505 

1755 

l 

25 

900 

460 

12650 

52 

325 

2080 

7475 

2 

36 

1656 

1012 

3107  0 

75 

600 

4600 

18975 

3 

46 

2530 

1771 

61985 

96 

920 

8096 

37191 

4 

55 

3465 

2695 

102410 

115 

1265 

1239? 

61985 

5 

63 

4410 

3724 

150822 

132 

1617 

17240 

92169 

6 

70 

5320 

4788 

203490 

147 

.1960 

22344 

1256Q5 

7 

76 

6156 

5014 

255816 

160 

2200 

27360 

159005 

0 

8.1 

6885 

6732 

302940 

171 

2565 

31977 

191862 

9 

83 

7480 

7460 

340340 

180 

2805 

35904 

216790 

10 

88 

7920 

8000 

364364 

167 

2992 

38896 

230238 

11 

90 

0190 

8281 

372645 

192 

3120 

40760 

240430 

12 

91 

8281 

0201 

364364 

195 

3165 

41405 

248430 

13 

91 

ai90 

8000 

340340 

196 

3165 

40768 

238238 

14 

90 

7920 

7480 

302940 

195 

3120 

38096 

21079C 

15 

03 

7480 

6732 

255816 

192 

2992 

35904 

191062 

16 

85 

6885 

5014 

203490 

187 

2805 

31977 

issues 

17 

01 

6156 

4788 

ISt  022 

180 

2565 

27360 

125605 

10 

5320 

3724 

102410 

i  7  l 

2280 

^  t  A  .1 

92169 

19 

70 

4410 

2695 

61985 

.160 

1960 

.17248 

61985 

20 

63 

3465 

1771 

31878 

147 

.163.7 

12397 

37191 

21 

55 

2530 

1012 

12650 

132 

1265 

8096 

10975 

22 

46 

1656 

460 

2990 

US 

920 

4600 

74  75 

23 

36 

900 

130 

96 

COO 

2080 

1755 

24 

25 

325 

75 

325 

585 

25 

13 

52 

1.17 

26 

| 

27 

27 

28 

.1630  118, 755 

101,790 

1 ,032,015 

3654 

47,502 

525, 915 1 

2,804,880 

lIn  computing  m(3)  in  Tabic  3,  the  weights  need  were  3  times  these  weight 
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TABLE  2  ( Continued) 


29  POINTS  30  POINTS 


i 

k  1 

PWj 

k=2 
pv  t 

k  ’3 

PWj 

k  =4 
pw, 

k  -1 

pWj 

k*2 

pw, 

k=3 

pw, 

k  =4 

pWj 

i 

.1 

14 

126 

819 

4095 

29 

203 

1827 

23751 

1 

2 

27 

351 

2925 

17550 

56 

567 

6552 

102375 

e 

3 

39 

650 

6500 

44850 

81 

1053 

14652 

2S32S0 

3 

4 

50 

1000 

11500 

88550 

104 

1625 

26000 

523250 

4 

5 

60 

1380 

17710 

148764 

125 

2250 

40250 

085500 

5 

6 

69 

1771 

24794 

223146 

144 

2698 

56672 

1538876 

6 

? 

77 

2156 

32340 

307230 

161 

3542 

74302 

1859550 

7 

a 

84 

2520 

39900 

335010 

176 

4158 

92400 

2413950 

0 

9 

90 

2850 

47025 

479655 

189 

4725 

109725 

2962575 

9 

.1.0 

95 

3135 

53295 

554268 

200 

5225 

125400 

3464175 

10 

11 

99 

3366 

58344 

612  6. 1.2 

209 

5643 

13856? 

3879076 

11 

12 

102 

3536 

61880 

6497  tu 

216 

5967 

148512 

4176900 

12 

13  . 

104 

3640 

63700 

662480 

221 

6188 

154700 

4331600 

13 

14 

105 

3675 

63700 

649740 

224 

6300 

156800 

4331600 

14 

15 

105 

3640 

61880 

612612 

225 

6300 

154700 

4176900 

15 

16 

104 

3536 

58344 

554268 

224 

6188 

148512 

3879876 

16 

1? 

102 

3366 

53295 

479655 

221 

5967 

130567 

3464175 

17 

in 

99' 

3135 

47025 

395010 

216 

''-a a  % 

v  v  *v 

n  Ciena  oc; 

Cm  \J  |  xj 

•» 

■  LO 

19 

95 

2050 

39900 

307230 

209 

5225 

109725 

2413950 

19 

20 

90 

2520 

32340 

223146 

200 

4723 

92400 

3859550 

20 

21 

84 

2156 

24794 

146764 

189 

4158 

74302 

1330076 

21 

22 

77 

1771 

17710 

88550 

176 

3542 

56672 

885500 

22 

23 

69 

1300 

11500 

44850 

16.1 

2090 

40250 

523250 

23 

24 

60 

1000 

6500 

17550 

144 

2250 

26000 

263250 

24 

25 

50 

650 

2925 

4095 

125 

1625 

14652 

102375 

25 

26 

39 

35.1 

819 

104 

1053 

6552 

23751 

26 

27 

27 

126 

81 

567 

1027 

27 

28 

14 

56 

203 

20 

29 

29 

29 

30 

Xw 

2030 

56,637 

041,464 

7,713,420 

4495 

100,680 

2,136,024 

52,451,256 

10 


TABLE  3 


BINOMIAL  COEFFICIENTS 


0 


1 

2 

3 

4 

5 

1 

2 

1 

3 

3 

1 

4 

6 

4 

1 

5 

10 

10 

5 

1 

6 

15 

20 

15 

6 

7 

21 

35 

35 

21 

8 

28 

56 

70 

56 

9 

36 

84 

126 

126 

10 

45 

120 

210 

252 

11 

55 

165 

330 

462 

12 

66 

220 

495 

792 

13 

78 

286 

715 

1287 

14 

91 

364 

1001 

2002 

15 

105 

455 

1365 

3003 

16 

120 

560 

1820 

4368 

17 

136 

680 

2380 

6188 

18 

153 

816 

3060 

8568 

19 

171 

969 

3876 

11628 

20 

190 

1140 

4845 

15504 

21 

210 

1330 

5985 

20349 

22 

231 

a  540 

7315 

26334 

23 

253 

1771 

8855 

33649 

24 

276 

2024 

10626 

42504 

25 

300 

2300 

12650  * 

53130 

26 

325 

2600 

14950 

65780 

27 

351 

2925 

17550 

80730 

28 

378 

3275 

20475 

98280 

29 

406 

3654 

23751 

118755 

30 

435 

4060 

27405 

142506 

31 

465 

4495 

31465 

169911 

32 

496 

4960 

35960 

201376 

33 

528 

5456 

40920 

237336 

34 

561 

5984 

46376 

278256 

35 

595 

6545 

52360 

324632 

36 

630 

7140 

58905 

376992 

37 

666 

7770 

66045 

435897 

38 

703 

8436 

73815 

501942 

39 

741 

9139 

82251 

575757 

Mil'll:' 


Given 


APPLICATION  OP  SMOOTHING  METHOD  TO 
SINE  CURVE  USING  UNEQUALLY  SPACED  DATA 


y  =  sin  x  , 

with  the  points  unequally  spaced  over  the  interval  0  <  x  <  1  ,  use 
of  the  divided  difference  method  described  in  the  text  yrelds 
results  tabulated  below.  First,  let 

y  =  fn  =  A0Pq  +  AjPjl  +  AsPa  +  . . .  +  AnPn  > 

where,  in  this  case  n  +  1  +  m  =  11,  n  =  4  and  Pk  =  xk  .  In  this 
example  u  =  x  . 

The  weights,  <*Ji  ,  are  chosen  by  the  following  modification  of  the 
weights  already  calculated  for  equally  spaced  data  containing  11 
points  (see  Table  2,  Appendix  A): 

<*>i  =  (ujt  -  u0)  wx 

wa  -  W+x  -  w2 


wi  =  ^uX+i+l  “  ul+l^  wi  > 


so  that  we  have 


1 


Then  j>roceeding  with  the  fit  we  obtain  the  results 


X 

y 

h 

(y  ■■  A4x4) 

(Pa  -  A3x3) 

Pi 

(Ps  -  A2x2) 

(P 1  “  AjX, 

0.00 

0.00000 

0  .00000 

0.00000 

0 .00000 

0.00000 

0.08 

0.07991 

0.07991 

0.08000 

0.07996 

0.00002 

0.21 

0.20846 

C  .20842 

0.21012 

0 .20986 

0.00001 

0.32 

0.31457 

0.31436 

0.32038 

0.31977 

-0.00001 

0.38 

0.37092 

0.37050 

0 . 38058 

0.37972 

-0.00001 

0.50 

0.47943 

0.47818 

0.50115 

0.49966 

0.00001 

0.61 

0.57287 

0.57011 

0.61182 

0.60960 

0.00003 

0.68 

0.62879 

0.62452 

0.68230 

0.67954 

0.00002 

0.82 

0.73115 

0.72213 

0.82344 

0.81943 

0.00000 

0.93 

0.80162 

0.78670 

0.93450 

0.92935 

0.00000 

1.00 

0.84147 

0.82152 

1.00527 

6.99931 

0.00001 

A4  =  0.01995 
A3  =  -0.18375 
A2  =  0.00596 
A1  =  0.99930 
A  =  0.00001 

O 

Therefore 

sin  x  ~  0.00001  +  0.99930  x  +  0.00596  x2  -  0.18375  x3 
+  0.01995  x4  ,  0  <  x  <  1  radian 


NOTE  ON  INTERPOLATION  WITH  TCHEBYCHEV  POLYNOMIALS 


When  Tchebychev  polynomials  are  used  in  equation  (1)  for  the  purpose 
of  truncating  to  a  lower  degree  the  u  should  be  chosen  so  as  to  moke 
the  balanced  oscillations  of  T(u)  occur  within  the  region  of  interest 
Thus,  if  accurate  results  are  desired  for  <  x  <  Xj  ,  u  should  be 
chosen  so  that 


u  =  -1  for  x  =  , 

and 

u  =  1  for  x  -  Xj 

This  is  particularly  important  in  the  process  of  inverse  interpolation. 
For  example,  suppose  one  is  given  the  following  table  of  values.  The 
problem  of  finding  the  solution  of  f(x)  =  0  will  obviously  involve  the 
investigation  of  a  neighborhood  near  x  =1.4. 


X 

<x) 

1.0 

-4.0000 

l  n 

±  .  a 

-2  .8224 

1.4 

-0.2464 

1.6 

4.1216 

1.8 

10.7136 

Linear  interpolation  using  the  points  (1.4,  -0,2464)  and  (1.6,  4.1216) 
gives  x  =  1.4.11  as  an  approximate  solution  of  f(x)  =  0.  Assume  that 
the  desired  solution  lies  between  x  -  1.40  and  x  =  1.42,  and  choose 
the  linear  transformation  connecting  x  and  u  so  that  u  =  -1  for 
x  -1.40  and  u  =  1  for  x  =1.42. 

The  following  table  can  then  be  used  for  determining  the  A’s  in  equa¬ 
tion  (1)  : 


X 

u 

T4(u) 

T3(u) 

Ta(u 

1.0 

-41 

22,592,641 

1.2 

-21 

1,552,321 

-36,981 

881 

1.4 

-1 

1 

-1 

1 

1.6 

19 

1 , 039,  681 

27,379 

721 

1.3 

39 

18,495,561 

237,159 

Using  these  values  of  Tchebychev  polynomials,  points  on  successively 
lower  degree  curves,  approximating  f(x),  can  be  obtained 


X 

f(x) 

^3 

’  f*3  **  A3T3 

fl  =  fa  -  AaTa 

1.0 

-4.0000 

1.2 

-2.8224 

-2.8224 

-2.73919275 

1.4 

-0.2464 

-0.2464 

-0.24639775 

-0.24753125 

1.6 

4.1216 

4 . 1216 

4.05999725 

3.24274375 

1.8 

10.7136 

10 . 7136 

A4  f  =  A4  A4  T4,  0.0384  =  30,720,000  A4 

A4  =1.25  x  10"9  f 

Because  of  the  magnitude  of  A4  the  term  A4T4(u)  is  neglected. 

A3  ffl  -  A3  A3  T3,  0.432  =  192,000A3 

► 

A3  =  C. 00000225 

A2  t3  =  A2  A2  Ta,  1.0136  -  1600A2 
A2  -  0.0011335 

Linear  interpolation  with  f1  gives  x  =  1.414184  as  an  approximate 
solution  of  f(x)  =  0.  This  can  be  used  as  a  first  approximation  for 
inverse  quadratic  interpolation  in  f o  ,  and  the  approximation  obtained 
from  f 2  can  be  used  for  starting  inverse  cubic  interpolation  in  f^. 

In  this  case 


f(x) 


x4 


+  3x3 


6x 


and 

r(x)  -  o 

for 

x  -  1  1.4142] 

so  that 


f, 

J. 

has  given  an  excellent  approximation  of  the  correct  solution. 


The  accuracy  of  the  various  approximations  for  f(x)  within  the 
interval  1.40  <  x  <  1.42  for  direct  interpolation  can  be  seen  from 
the  A's: 


If  "  f3l  <  |Aj  =  1.25 
If  -  fsl  <  |Aj  +  I Ay | 
if  -  fj  <  Ia41  +  |a3| 


10“  9 

1.25  x  10" 9  i-  2.25  x  10” 6 
lAs|  -0.00113575 


